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ABSTRACT 


(•OOc. 


_\J  V  -4  v  -2 

Electrical  conductance  data  in  the  concentration  range  10  to  10 

for  LiBF^in  1,2-DME  are  interpreted  by  the  Fuoss-Kraus  theory  yielding 

values  of  the  ion-pairs  and  triple-ions  formation  constants. 

Complex  dielectric  permittivities  for  LiBF^lin  1,2-DME  at. 25*^0  in 
the  frequency  range  0.9-81  GHz  are  reported.  The  data  are  interpreted 
by  two  Debye  relaxation  processes,  one  for  the  solvent  and  one  for  the  solute. 
The  latter  is  interpreted  as  due  to  the  diffusional  rotation  of  ion-pairs. 

Ultrasonic  absorption  data  for  LiBF^in  1,2-DME  at  25^C  in  the  con¬ 
centration  range  0.28-0.55M  and  frequency  range  10-550  MHz  are  reported. 
The  data  are  described  by  a  single  Debye  relaxation  process  and  are  in¬ 
terpreted  as  due  to  an  ion-pair  dimerization  equilibrium  leading  to  forma¬ 
tion  of  quadrupoles.  Theoretical  expression  (similar  to  the  ones  of  Fuoss 
and  Bjerrum  for  ion-pairs)  are  developed  for  the  dimerization  of  ion  pairs 
to  quadrupoles.  To  this  end,  a  point  dipole -point  dipole  potential  has  been  used. 
The  dependence  of  K'p  and 'i6^?from  tfte  dipole-dipole  separation  distance 
is  illustrated.  Further,  the  thermodynamic  increments  ASq,  AHq 

and  AVg  have  been  calculated^  The  latter  quantity  is  compared  to  the 
corresponding  parameter  foun<^  during  the  present  and  previous  ultrasonic 
research.  \  '  . 

^  W<v  V  -  ,  V.  3  ••  ^  *  1  T  ^  •  ,  ■.  c  ,L  \  v.  , 


Itu,  ^ 


*The  ultrasonic  data  for  LiBF4  in  1,  2-DME  have  been  already  given  in 
U.  S.  Army  Report  No.  3,  dated  June  30,  1983,  under  Grant  No.  DAAG- 
29-82-K0  048.  They  are  regiven  here  for  the  sake  of  logic  of  the  present 
report  that  will  be  basically  a  paper  to  be  submitted  for  publication. 


INTRODUCTION 


1,2  Dimethoxyethane  (1,2-DME)  is  known  to  chelate  lithium  as  well  as 
•odium  ions.  LiC tO^  [  1  ]  and  ULAsF^  [21  have  been  investigated  in  this 
solvent  by  ultrasonic  and  dielectric  relaxation  and  conductance  methods. 
LiBF^,  another  electrolyte  relevant  to  a  battery  construction,  has  now 
been  investigated  by  the  same  methods  in  1,  2-DME  in  order  to  establish 
whether  a  general  trend  of  association  exists  in  different  solvents  (a 
sequence  has  already  been  found  for  the  three  above  electrolytes  in 
2-Methyltetrahydrofuran  [3]). 

Enough  systems  have  shown  formation  of  quadrupoles  in  ethereal 
solvents  to  justify  interest  in  theoretical  expressions  aimed  at  calculating 
thermodynamic  parameters  relevant  to  the  ultrasonic  research.  To  this 
end,  using  a  point  dipole-point  dipole  potential,  theoretical  expressions  of 
quadrupole  formation  constants  and  thermodynamic  increments  as  ASq, 

AH  ,  AV™,  and  AV  have  been  obtained  as  shown  below.  In  some  cases 
already  investigated,  a  process  of  dimerization  consisting  of  two  steps 
appeared  to  exist. [3].  For  these  cases,  the  second  step  of  dimerization 
is  interpreted  theoretically  by  the  Kirkwood  expression  of  dipole  associa¬ 
tion  [4].  - 

EXPERIMENTAL 

The  equipment  and  procedure  for  the  ultrasonic  dielectric  and  conduc¬ 
tance  work  have  been  described  elsewhere  [31.  LiBF^  (Aldrich)  was  re¬ 
dried  in  vacuo  (-  1  torr)  at  -  70°C  overnight.  1,2-DME  has  been  purified 
by  distillation  at  reduced  pressure  over  sodium  and  benzophenone.  Solu¬ 
tions  were  prepared  just  before  each  run.  Runs  performed  with-solutions 
aged  3-4  days  gave  results  which  were  outside  the  experimental  error 
by  a  systematic  amount.  These  results  were  discarded. 
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RESULTS 


Fig.  la  reports  the  equivalent  conductance  A  vs.  the  concentration 
c  in  a  double  log  plot  for  LiBF^  in  1, 2-DME  at  25°C.  In  the  same  plot, 
previous  results  [2]  for  LiAsF^  in  the  same  solvent  show  that  LiAsF^  is 
a  stronger  electrolyte  in  this  solvent.  This  confirms  the  same  sequence 
found  in  2-Methyltetrahydrofuran  [3]. 


Fig.  2a  reports  the  coefficients  of  the  real  part  and  of  the  imaginary 
part  e#  and  e*  of  the  complex  permittivity  e  =  plotted  vs.  the 

frequency f  for  a  representative  concentration  of  LiBF^  in  1, 2-DME  at 
25°C.  The  solid  lines  are  the  calculated  Debye  functions  for  two  relaxa¬ 
tion  processes: 


€o"€».  Sao  "eoo 
£  '  =  ?  +  —  - r—  +  C 

1  +(f/fj)  l+(f/f2) 


oc. 


(1) 


c'-e'  =  (e_-e„  ) 


f/fl 


+  (e«  ) 


f/f. 


'  /  c  '  “  eo  ao  '  . 

i  i+(f/fxr  i  2  i+(f/f2> 


(2) 


with  parameters  e  .  G  ,  €  ,  f , ,  f  •  e*  is  the  conductance  contribution 

o  j  l  w  x 

1  8  x 10* 

to  the  loss  “ — £ - -  ;  x  is  the  specific  conductivity  of  the  solution. 

Table  I  reports  all  the  dielectric  parameters  and  x  for  the  concentra¬ 
tions  of  LiBF^  investigated  at  25° C. 


Fig.  3  reports  a  representative  plot  of  the  quantity  ct/x  vs.  the  fre¬ 
quency  f  for  LiBF  .  in  1,  2-DME  at  25°C.  ais  the  sound  absorption  coefficient 
’  _  1 

expressed  in  neper  cm  .  The  solid  line  is  the  fitted  Debye  function  for 
a  single  relaxation  process: 


JL  = 


A 

(1  +(f/£r)2 


+  B 


with  A,  B  and  L  the  relaxation  parameters.  In  particular,  A  (=  2umax/fru) 
contains  all  the  thermodynamic  information  relevant  to  the  observed  proces 
^max  *8  maximum  absorption  coefficient  per  wavelength  at  f  =  fy,  where 


£  is  the  relaxation  frequency.  B  is  the  background  absorption  at  f  »f  ; 

*  2  * 
namely,  the  value  of  n/t  at  high  frequency  with  respect  to  where 

the  degrees  of  freedom  of  the  observed  process  have  been  frozen  (relaxed). 

u  is  the  sound  velocity.  Table  II  reports  the  calculated  parameters  for 

the  concentrations  of  LiBF^  in  1,2-DME  investigated  at  25°C. 

DISCUSSION 

a.  LiBF4  in  1,2-DME 

The  conductance  data  have  been  treated  by  the  Fuoss- Kraus  theory 
of  triple- ions  formation  [5  ] 

_  K  A?*K_  , 

Ag(c)v^=Tz:  +  “zr(i-f  )c  ,  (4) 

,  -/K,  VKp  Ao 

where  g(c)  is  a  factor  that  lumps  together  all  the  interionic  interaction 

T  +  - 

terms  Aq  and  Aq  are  the  limiting  equivalent  conductance  of  Li  +  BF^ 

and  of  the  two  possible  triple  ions  (L^BF^  +  Li(BF4)2),  respectively. 

Kp  and  K._,  are  the  ion  pair  and  triple  ions  formation  constants,  respectively. 

The  symmetrical  approximation  of  considering  the  two  possible  formation 

constants  of  triple  ions  equal  to  each  other  has  been  retained. [  5  ] . 

Fig.  lb  reports  the  plot  of  Ag(c)*/c  vs.  c(l  according  to  Eq.  (4). 

Ao 

The  solid  line  has  been  calculated  by  linear  regression.  It  gives  the  de- 

2  7 

termination  coefficient  r  =0.989,  intercept  =  3. 921  x  10  ,  slope  =  1. 3122. 

From  a  previous  work  [3^]  the  value  of  the  Walden  product  A^ri  =  0.564  for 

LiBF.  was  calculated.  Hence,  for  this  electrolyte  in  1,2-DME  one  estimates 
*  - 1  2-1 

A0  =  136.6D  cm  eq  .  From  the  above  results  and  the  arbitrary  position 

T  2 

AQ  =  "j  AQ  (consistent  with  previous  work)  [  1 , 2,  3 "]  one  calculates 

K  =  1. 2  x  10?M-  1  and  K„=50M"1 
p  T 

It  is  interesting  to  compare  these  figures  with  the  corresponding  ones  [21 
for  ULAsF^  in  1, 2-DME  K  -  1  x  10^  M  ^  and  K^,  =  28  M  * .  One  may  con¬ 
clude  that  the  reason  that  LiAsF^  is  a  stronger  electrolyte  in  1,2-DME  with 


4- 


respect  to  LiBF^  is  due  to  the  difference  of  two  orders  of  magnitude  in 

the  value  of  K  rather  than  to  differences  in  the  triple  ion  formation 
P 

constants,  K^,.  It  is  also  relevant  to  compare  the  conductance  results 

of  the  two  electrolytes  LiAsF,  and  LiBF.  in  the  two  solvents  1,2-DME  and 
$  o  *±  7—1 

2MTHF.  For  the  latter  solvent  [3^»  LiAsF,  showed  K  =3.3x10  M  , 

-1  ®  P 

K^,  =  39M  .  Again,  the  gain  in  orders  of  magnitude  of  the  equivalent 

conductances  of  both  electrolytes,  when  dissolved  in  1,2-DME  with  respect 

to  2MTHF,  seems  due  to  the  differences  in  the  values  of  K  rather  than 

P  + 

the  K.j.'s.  On  a  molecular  scale  it  is  likely  that  chelation  of  Li  by  1,2-DME 
shifts  equilibrium  (9)'  (see  below)  toward  the  left  thus  decreasing  Kp. 

On  the  other  hand,  the  steric  hinderance  caused  by  the  presence  of  a  -CH^ 
group  adjacent  to  the  ethereal  group  of  the  2MTHF  may  cause  solvation 
hinderance  which  favors  the  anion  as  a  competitor  for  the  first  coordination 
shell  around  Li+.  This  may  shift  equilibrium  (9)  toward  the  right  increas¬ 
ing  Kp.  It  is  also  important  to  realize  that  the  major  proportion  of  the 
electrolyte  exists  as  ion- pairs  with  only  a  minor  proportion  as  triple  ions 
(neglecting  quadrupoles  which  are  revealed  by  ultrasonic  relaxation  as 
shown  below).  On  the  basis  of  the  conductance  information  alone  for  the 
equilibrium  [  5  1 ,  A  +  Bf^,  AB  +  B ^  AB2,  (AB)  =  (l-a~3a^.)  c«c,  (AB^)  Ssa-j-c 
(B)  s*  cl  c  ,  a«(VKpy^  r1,  aT  =  (KTA/iT)cl/2,  hence  (AB^  (KTA/#)c3/2. 
With  the  above  figures  K  =  1. 2  x  10^  M“  *  and  K^,  =  50  M  atc  =  0.1  , 

(AbZ)  =  4. 55  x  10-4M  or(AB~  /c)%  =  0.46. 

r  2  60+1--, 

Fig.  4  depicts  the  quantity  I  ( c  ^ -  I  vs*  c  according  to 

1  “o 

the  Bottcher  function  [6  3: 


(e-e  )..«ttl«cx10 

0  1  (1-af) 


3k  T  2e  +1  * 

o 


where  the  polarizability  a,  reactor  field  factor  f  term  (1-af)  generally 
of  the  order  of  0.9  has  been  equated  to  1.  In  the  above  4  is  the  apparent 
dipole  moment  of  the  rotating  polar  entity  which  increases  the  polarization 
of  the  solution  and  makes  e  >  .  The  solid  line  in  Fig.  4  has  been 

calculated  by  linear  regression  (5u%  stat.  weight  to  the  origin)  giving 

r2  » 0. 98,  int  =  0.031  4  0.  16,  slope  =  15,  255  ±  2743  from  which  4  =  (15.8  ±1.4) 

"18 

x  10  e.s.u.  cm.  By  taking  a  rigid  dipole  moment  4  =  ae  where  e  is  the 


2MTHF  =  2-Methyltetrahydrofurao . 
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electron  charge  a  -  3.3  ±  0.3  x  10  cm.  This  figure,  when  com¬ 

pared  to  the  one  obtained  for  LiAsF,  in  the  same  solvent  from  a  previous 

-8  ® 

work  [  2  3  a  =  (3. 8  ±  0.  04)  x  10  cm,  suggests  a  smaller  charge  to  charge 

separation  of  the  rotating  entity  for  LiBF^  with  respect  to  LiAsF^.  This 

may  simply  reflect  the  molecular  dimensions  of  the  two  anions.  On  the 

other  hand,  LiCfO  .  in  1, 2-DME  gave  [13  pi  =  14.  7  x  10  esu  cm  and 
_  s  J* 

as  3. 1x10  cm.  The  following  sequence  for  the  apparent  dipole  moments 
and  charge  separation  seems  to  appear  LiAsF^  >  LiBF^>  LiC2  O^.  As  the 
crystallographic  radius  of  Cl  is  larger  than  BF^ ,  it  appears  that  other 
factors  i.e.,  solvation  and  the  different  ability  of  the  anions  to  penetrate 
the  solvation  shell  of  Li  .play  a  role  in  the  found  sequence  of  the  a^'s. 

It  is  also  noteworthy  that  for  all  the  above  electrolytes  in  1,  2-DME  fr  t  2.  0 
GHz,  hence  the  decay  time  of  the  polarization  is  constant.  Neglecting  • 
differences  between  this  latter  quantity  and  the  molecular  relaxation  time 
f  and  retaining  the  Debye  expression  [73  for  a  spherical  dipole 


dtra" 
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f  a  — ,  ti,  it  results  that  a  =4.0x10  cm  expressing  the  average  solvated 
Kl  T 

radius  of  the  rotating  entity. 


The  ultrasonic  data  will  now  be  discussed.  Linear  dependence  of  fr 
on  c  and  nonlinear  concave -up  dependence  of  on  c  suggest  [13  a 

dimerization  process  of  the  type; 


with  A  symbolizing  the  LiBF^  ion-pairs  and  B  the  dimers  (LiBF^)^.  The 
following  expressions  are  valid  for  this  process  [13: 


T*1  =  2Trfr  =  4k£(A)  +  kr  ,  (6) 

(AV  )2  K  (A)2 

*m=3r8-R%“  •  (?) 

q 

where  k,  and  k  are  the  rate  constants,  AV  is  the  isoentropic  volume 

1  r  2  - 1  3 
change  for  the  process,  g  =(pu  )  is  the  adiabatic  compressibility  of 

8 
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the  liquid,  p  the  density  andK^  =  k^/k^  the  equilibrium  constant  of  de- 
merization. 

To  a  first  approximation  if  is  not  large  (-1  -  3M  *)  one  may 

approximate  (A) «  c  taken  as  the  concentration  of  ion  pairs,  the  major 

species  present  in  a  solvent  of  this  permittivity  (e  «  7).  A  plot  of  t  * 

vs.  c  is  depicted  in  Fig.  5a.  Linear  regressions  give  slope  =  9.68x10°, 

8  2 

intercept -2. 95  x  10  ,  determination  coefficient  r  =0.84  (r  is  the  cor¬ 
relation  coefficient).  Then,  k,  =  2. 4x10®  M  k  =  2.95xl0®S~*  and 

-1  1  r 
K  =0.8  M  .  Equation  (7)  can  be  transformed  to  the  Farber's  function 


-*8 


1 13* 


.ns. 


2fi  4  RT 


VAVs> 


(AVs) 


*8) 


-2 

The  Farber  plot  is  shown  in  Fig.  5b.  It  gives:  slope  =  8.65x10'  , 

-5  2 

intercept  =  0.90x10  ,  r  =0.81  from  which  one  calculates  AV  =6.8 

3  —  1  ® 

cm  /mol  and  =  2.4M  ,  of  the  same  order  of  magnitude  of  the 

figure  K  =  0.8  M  1  obtained  above.  Fig.  5c  and  Table  II  report  the 

9 

calculated  and  the  experimental  4  's  testing  the  internal  consistency 


-1 


of  the  calculated  parameters  AV  and  K  =  2.4M 

s  q 

Comparing  the  above  ultrasonic  results  with  the  ones  obtained 
for  LiAsF^  and  LiCfO^in  DME,  one  has  to  conclude  that  ultrasonic 
relaxation  seems  able  to  detect  one  of  the  following  equilibria,  depend¬ 
ing  on  the  system: 


k 


U+  +  B' 


Li+. 


-o 


B  LiB 
k-l 


(9) 


2 LiB  ;=±  (LiB). 


B  =  ligand 


namely,  one  of  the  steps  leading  to  the  contact  ion  pair  LiB  or  to  the 
dimer  (or  quadrupole  (LiB^h  Generally,  in  media  of  low  permittivity 
the  first  step  of  the  scheme  above  is  not  detectable  because  of  the  small  concentra¬ 
tion  of  the  free  ones.  It  appears  that  for  LiAsF^  in  1,2-DME  the  step 


7- 


.  •  -  ■  *  '  •  -k  . -  .  -  . 


-V  'V  ' 


Li*...B  LiB  was  detected  [2 3,  whereas  for  both  LiBF^  and  LiC£ 

[13  only  the  dimerization  is  visible.  This  could  reflect  the  extent  of 
association  of  the  different  electrolytes,  LiAsF^  being  associated  only  to 
the  form  LiB,  dimers  being  in  a  minor  quantity.  Similarly  in  a  poorer 
solvent  as  2-Methyltetrahydrofurane,  LiAsF^  shows  formation  of  quadrupoles 
by  a  single  step  [3],  LiBF^by  a  double  step  [3^  and  LiCfO^  only  by  a 
first  order  process  [8]  according  to  the  sequence 


2  LiB  ^  (LiB...  LiB) 
(LiB...  LiB)  ~±  (LiB) 


(10) 


This  would  show  that  LAC2  is  mainly  in  the  form  of  quadrupoles  of  the 
outer  sphere  to  inner  sphere  type  [83  and  the  sequence  of  association 
•LiAsF^  LiBF^<  LiCf  is  confirmed.  At  this  point  of  the  research  it 
was  felt  that  enough  systems  had  shown  presence  of  quadrupoles  and 
enough  parameters  as  and  AVg  had  been  determined  to  justify  interest 
in  calculating  the  order  of  magnitude  of  these  quantities  on  a  theoretical 
ground.  This  shall  be  attempted  in  the  sections  to  follow. 

b;  Calculation  of  the  Quadrupole  Formation  Constant  for  Dipole- 
Dipole  Interaction 

We  shall  start  our  calculation  with  a  Fuoss-like  expression  for 
dipole- dipole  interaction.  As  already  shown  by  Saar  et.  al.  [  9  3,  one 
can  paraphrase  the  Fuoss  expression  for  the  association  constant  of 
simple  ions  [5  3  to  a  dipole-ion  pair  to  the  one  of  dipoles  to  a  quadrupole 
by  writing 

KF  =  "3oorexp('rr)  f11) 


Where  r  is  now  the  separation  distance  between  the  centers  of  the  two 

2  3 

dipoles  of  moment  u  and  U  is  the  dipole -dipole  potential  U  =  -|i  /er  , 
with  e  the  permittivity.  One  then  writes 

K  =  lniill  exp(  )  =  K  eb 


1T1T1 


(12) 


with  Kq  =  4TrLr  /3000  and  b  =  (j  /er  kT. 

Refinements  to  this  equation  could  be  added  as  done  by  Jagodzinski  [10] 
for  the  association  of  ions  to  ion-pairs.  One  can  consider  the 
change  in  the  number  of  degrees  of  freedom  in  going  from  two  dipoles  to 
a  quadrupole.  By  using  the  classical  equipartition  principle  (and  reasoning) 
that  each  ion  pair  will  have  3- translational,  2- rotational,  and  1-vibrational 
degree  of  freedom)  one  writes  for  the  non- electrostatic  internal  energy 
[10]  of  the  two  ion  pairs 

2Unei  =  2[|rT  +  RT  +  Rt]  =  -y  RT 

$ 

The  quadrupole,  on  the  other  hand,  taken  as  a  nonlinear  molecule,  namely, 
a"  ±  t  "  entity,  will  have  3-translational,  3-rotational,  an  (3N-6)  =6 
vibrational  degrees  of  freedom.  Then 

u  .  =-|rt+-|rt+  6rt  =  —  RT 
net  2  2  2 


The  net  result  of  the  formation  of  a  quadrupole  from  two  dipoles  will  be  a 
change  in  non- electrostatic  internal  energy  per  mole  of  AUne  =  2RT, 
whereas  the  total  energy  change  per  mole  will  be 


AU  =  AU  +  U  =  2RT  -  L 
ne  e*  „  3 

S  T 

with  JL  the  Avogadro's  number.  Then 

=  45feLe-2e^'/er3kT  =  K*eb 

o 


K 


3000 


(13) 


(14) 


with 


4TrT_.r^  -2  ,  ,  V? 

K  =  — e  and  b  =  — i 


3000 


3.  _ 
cr  kT 


*v* 

Although  ultrasonic  relaxation  has  detected  quadrupoles,  dielectric 
relaxation  for  Li  salts  in  the  same  solution  has  never  detected 
these  species.  Hence  they  must  have  a  zero  or  negligible  dipole  moment. 


$ 

Table  III  reports  the  values  of  K_  and  K  vs.  the  interaction  distance 

F  _  jg  F 

r  for  two  dipoles  of  moment  (a  =  14.7  x  10  esu  cm  immersed  in  a  medium 

of  permittivity  e  =  7  at  T=  298. 2K.  It  may  be  seen  that  both  functions  go 

through  a  minimum  by  increasing  r.  This  is  self  evident  by  observing  that 
*  3  l 

K_  or  K  depends  on  the  product  r  e1' r  ,  namely,  the  product  of  the  two 
it  2 

functions,  one  increasing  and  one  decreasing  with  r  .  The  position  of  the 

minimum  can  be  found  by  maximizing  Eqs.  12  or  14  .  Then  in  the  case  of 

Eq.  14: 

kT]=° 

$ 

which  gives  the  distance  at  the  minimum  of 

/  u2  N1/3  -8 

rm  ~  VeltT )  =9.1x10  cm  (15) 

It  may  be  useful  to  recall  [  1 1 ]  that  Kq  is  a  measure  of  the  excluded  volume 
due  to  the  presence  of  neutral  dipoles  or,  in  other  words,  the  probability  of 
association  increasing  with  the  volume  occupied  by  neutral  dipoles. 

Some  researchers  [12]  in  the  past  have  expressed  preference  for 
the  Bjerrum  theory  of  ion-pair  association  over  the  Fuoss  one  based  on  the 
fact  that  the  former  seems  to  follow  the  data  for  the  formation  constant  with 
€  more  closely  than  the  latter  one.  It  appears  natural  then,  at  this  point, 
to  develop  a  Bjerrum-like  theory  for  quadrupole  formation  and  compare  it 
with  Eqs.  12  and  14. 


Following  Bjerrum  [13],  the  probability  of  finding  a  dipole  at  a 
distance  r  from  a  given  one  will  be  given  by  the  radial  distribution  function: 


dP  = 


Lc 

1000 


4rrr2  dr  e'U/kT 


(16) 


where  Lc/1000  is  the  concentration  of  dipoles  expressed  as  molecule  cm  , 
2 

4rrr  dr  is  the  spherical  volume  element  of  radius  r  and  thickness  dr,  U  is 

2  3 

the  dipole- dipole  potential  U  =  - p  /er  .  Then, 


-10- 


The  probability  function  will  have  a  minimum  with  r  as  it  depends  on  the 
1  /  2 

product  e  x  r  .  Therefore 


dP 

dr 


Ar 


Lc 

1000 


4n-^[(exp 


u 


er3kT 


which  leads  to  the  condition 


q  = 


m 


_f3  Jilll/3 

-  L2  ekTj 


(18) 


Notice  that  the  minimum  in  the  distribution  function  corresponds  to  the 
position  |j  /eq  =  2/3  kT,  namely,  to  the  ratio  of  the  dipole-dipole  energy 
to  kT  equal  to  2/3.  Following  Bjerrum  [13],  we  shall  limit  the  definition 
of  associated  dipoles  to  the  distances  r  in  the  range  from  a  to  q,  with  a 
the  minimum  approach  distance  of  two  dipoles  of  moment  4.  The  probability 
that  two  dipoles  are  associated,  namely,  the  degree  of  association  1-c  will 
be  given  by 

1-0  ■  /V2/CT3kT)r2dr 

a 


Call  the  quantity  Y  =  |j2/er3kT  and  b  =  [j^/ea^kT.  Then,  3r^dr  =  -( p^/ekT) Y  2dY 

and 


1-0 


2 

4ttLc  u 
1000  3ekT 


since  for  r  =  a,  Y  =  b  and  for  r  =  q,  Y  =  2/3.  Then  the  equilibrium  constant 
expressed,  in  dilute  solutions,  after  Bjerrum  [13],  as  Kg.  =  1-o/c,  will 
be  given  by 


4ttL  UZ 
1000  l3ekT 


Y 

e 


Y_2dY 


b 


(19) 


or 


K„.  = 


_  4rrLa' 


Bj  3000 


b/ 

2/3 


Y  -2 
e  Y  dY 


(20) 


r  Y  -2 

The  integral  Q  =  J  e  Y  dY  can  be  solved  by  expanding  in  a  convergent 

2/3  234 

Y  Y  Y4*  YJ  Y 

series  e  ,  namely  e  =  1+Y  +  ’3*“  +  “41“  •••  •  Then, 


«»■/  [Vv  +ii  +  tr  +  ir  •••]dY 

0.67  Y 


or 


Eq.  21  can  be  written  in  a  more  compact  way  as 


.+ 


1 


(n+1).  n 


r  —  Yn"1 


0.67 

(21) 


Q  =  1.8930 +  lnb  -  £  +  t  7^+TjTn  [^”-(0.67)“] 


(22) 


It  appears  that  in  a  medium  of  permittivity  e=  7,  for  dipoles  of  moment 
ji  =  14. 7x10  esu  cm  at  T  =  298 . 2  K,  a  sufficiently  approximate  value  of 

_g 

Q  can  be  obtained  for  a  >  5  xlO  cm  by  limiting  n  in  Eq.  12  to  n=10. 

Table  III  gives  Kgj  as  a  function  of  r  =  a  and  Fig.  6  compares  Eq.  20 
with  Eqs.  12  and  14  (although  the  comparison  with  Eq.  14  is  not  con¬ 
sistent,  the  evaluation  of  AU  ,  having  been  omitted  in  the  evaluation 

nonet 

of  Kd.).  Equation  20  predicts  obviously  K-..  =  0  for  b  =  0.67,  which  cor  res - 
tij 

ponds  to  a  =  q.  At  this  point  the  same  refinement  to  K^.,  namely,  the  intro¬ 
duction  of  a  non- electrostatic  potential  change  can  be  applied  to 

Kt»*.  The  matter  reduces  to  write  the  total  potential  change  per  molecule 
of  dimer  2 

U  =  2kT  - 

er 


12- 


»-«•& W'  S'  7"  kT)^ 


which  becomes 


4ttLc  -2 
3000  e 


b  /  eY  Y-2  dY 
2/3 


The  above  expression  can  be  written  in  the  form 

KBj  =  Kob« 

£ 

Table  III  gives  the  values  of  Kgj  for  various  values  of  r.  The  same  quantities 
are  depicted  in  Fig.  6. 

**  Calculation  of  the  Thermodynamic  Parameters 


We  shall  now  proceed  to  calculate  the  thermodynamic  parameters 
relevant  to  the  ultrasonic  work  namely  the  isoentropic  AVg  °*  tlle  chemical 
process  of  ion-pair  dimerization  to  form  quadrupoles.  As  the  AVg  is  related 
to  the  AV^,  the  isothermal  volume  change  of  the  process  by  the  relation 


AVs  =  avt  "  7c~  AH  <24> 

p 

it  vnll  be  necessary  to  calculate  AH°  as  well.  In  the  above  0  =dinV/c>T  is 

the  expansivity,  p  and  c  the  density  and  specific  heat  of  the  liquid,  re- 

P  m 

spectively.  .  • 

The  calculation  for  AV^  will  follow  analogous  lines  of  the  estimation 
of  AV.J,  for  the  formation  of  ion-pairs  from  free  ions  after  Hemmes  [  14]  . 
The  expression  for  the  formation  constant  of  quadrupoles  to  be  bound  to  the 
AV.J,  has  to  be  expressed  in  molality 

,dlnK _ .  AVm 


-Vs)  ■- 


and  Km  is  related  to  the  equilibrium  constant  expressed  in  molarities  by 
the  relation 


Therefore 


in  K _  =  InK  +  in  p 

m  r 


/  dinKN 
V  dp  L-- 


fdln  o 
^  dp 


) 


KT  "  $  t 


with 


/ din  o  \  /  dinVN  _  1  /jWA 

V  dp  )-  V  dp  ;  "*P  T  'vlspl 


where  is  the  isothermal  compressibility  of  the  liquid.  We  shall  use 
for  K,  Eq.  14 

cr  kT 

,dinK*  2  ,  -1  .  2  ,  , 

/ _ F\  _  u  de  de  _  u  .  -2  de 

^  dp  L~  At  de  dp  'At  dp 


,dinK, 


dp  T  r3GkT  V  dp  T 


/dine  N 

V  dp  Jr 


Let  us  call  after  Fuoss  r  =  a  as  the  contact  distance  of  collision  between 
two  dipoles  when  they  are  defined  as  associated.  Then  by  calling 


3  i  ^ 

a  ekT 


•  fee; 


In  order  to  have  an  evaluation  of  AV^,  m  1,2-DME  the  following 
approximate  calculation  will  be  performed.  At  25°C  for  1,2-DME 

the  density  p  =  0.8614  gr/cm-*  [15  3  an(*  the  sound  velocity  u  =  1 .  18  x  10^  cm/ 

2-1  -12-1  2 
•  ec.  The  isoentropic  compressibility  8g  =  (pu  )  =83.37x10  gr  cmsec  . 

We  shall  approximate  0  0-,*  Also,  for  toluene  Fuoss  et  al.  [  16]  reported 

C  (1  at  m)  =  2.375  and  e  J  (5  80  8  at  m)  =  2.  775,  hence  (df  ne/dp)-_  = 


2. 7  xlO  ^  at  m  * , 


By  retaining  this  value  for  1,2-DME,  the  figure  of  Table  IV 


for  AV,j,  as  a  function  of  the  distance  a^  have  been  collected. 

We  shall  now  proceed  to  the  calculation  of  AH°.  From  the  expression 
*  b 

K  =K  e 
F  o 

AG  =-RTinK  =-RTfnK*-  RTb 
o  F  o 


*dAG  *  .  ^ 

ASo  =  -  (if)  *  R1”Ko  tRbtRTM 


■2>  -jiii  ,e-iT-i,.jiirT-1»£:^.  +  e-1ieL_i 

3T  3,  3T  e  1  "  3,  56  5T  6  5T  J 


3k  =Jiir.i€-2  le  _±t-2-i 
5T  .3.  L  T  5T  €  J 


'Sine  J_ 

-v  m  i  m 


a  ckT 


3tl  = 
aT 


X.  rafne  .  1  1 

•blTaT  +tj  • 


Then, 


AS  =  RfnK*  +  Rb  -  RTbFM^r  +41 


AH  =  AG  +TAS  =  -RTfnK*-  RTb  +  RT  inK*  -  RT2b 

O  o  o  o  o  ST 


*»„  *  -RTb[l  +  ^r] 


In  Table  IV  the  calculated  AH  's  are  reported.  The  parameters  used 
were  e  =  7,  4  =  14.7x10  esucm,  while  the  value  of  dfne/dT  =  5x10 
has  been  taken  from  the  data  of  Jagodzinski  [  10  3  for  THF,  assuming  a 
similar  temperature  depence  of  e  fori,  2-DME.  This  leads  to  (d£ne/d£nT) 
1.491.  The  above  calculation  is  admittedly  approximate,  but  as  we  are 
looking  for  an  order  of  magnitude  of  AVg  it  appears  adequate  for  the 
purpose. 

The  calculation  proceeds  with  the  estimation  of  0,  the  expansivity. 
From  the  densities  of  THF  with  temperature  of  Jagodzinski  [10],  one 


obtains 


o  _  .9 V  _  Sf  nV  _  Sf up  ,  in  ^ 

0- V  ST"  ST  "  ST  1*18xl°  K 


between  25  and  -15CC.  The  specific  heat  of  diethyl-ether  [17]  Cp(30°Q 
2.30  Joule  /K,  gr.  Then  the  quantity 


-  kll  H°T(T  =  5.817x10 
pc  0.88x2.30 


-4  cm 
Joule 


5. 817x1 0”^  cm8/ergs 


The  term  (0/ pc  )  AHq  expressed  in  cm  /mole  is  then  reported  in  Table  IV 
together  with  the  valculated  AV  's  according  to  Eq.  27. 

So  far  for  LiBF^  in  1 , 2-DME  and  for  LiCIO^  in  1 , 2-DME  the  values  found 

for  the  AV  's  were  of  the  order  of  7-10  cm^/mole.  Judging  from  Table  IV,  one 

would  estimate  a  collision  distance  between  the  two  dipoles  of  the  order 
-8 

of  6-7  x  10  cm,  a  very  reasonable  range  of  figures  for  any  model  envisaging 

-  8 

two  dipoles  of  charge  separation  2.5  -  3x  10  colliding  to  form  a  quadru- 
pole  .  In  addition,  the  experimental  value  of  the  quadrupole  formation 
constant  has  been  found  to  be  in  the  range  of  1-3  M  Judging  from  Table 
111,  equating  =  Kp  would  again  lead  toward  a  =  6-7x10  0  cm. 


- 16- 


It  appears  proper  at  this  point  to  repeat  the  above  sequence  of 
calculations  with  the  Bjerrum-like  function  20  of  quadrupole  formation 
constant. 

3  2^ 

K  -  v  kn  -  4frLa  u  r  Y  -2 , 

KBj  '  Kob  Q  -  3000  3,  _  J  e  Y  dY 

“  kT  2/3 


Given  the  relation 


dinK  A  V 

_ _ Bi  _  _ T 


RT  "  PT 


Also 


ana 


lnKBj.  =  fnKQ  +  i  nb  +  f  nQ 

_  dlnb  db  dfnQ  dQ  db 
dp  db  dp  +  dQ  db  dp 


db  _  _  u2  de"1  ds  _  .  dine 
dp  a3kT  de‘  dp  dP 


Therefore 


d/„K 


dlnKD. 
_ Bi 


Ti  +1  eb  1  din  e 

L  Q  b  J  dp 


AVX**T{[1+i%-]^-BT} 


As  done  above,  in  order  to  calculate  AV  one  needs  an  estimation  of  AH 
Then:  °* 


AG  =-  RT  in  K_  . 

o  XJJ 


AS  = 
o 


f  =  RinK  +  RT^f  +  ^  5b  ]£ 


as  found  above 


db  uTI  +  dfne"l 
dT  "  "  LT  +  dT  J 

AS0  =  R«nKB.  -  RTb[I  +  £ 

ASq  =  RinKfij  -  R[l  +  q  ~]  [l  +  d£nT  ] 


and  AH  =  AG  +  TAS  gives 
o  o  o 

k][‘ 


dine 

dfnT 


Eqs.30  and  31  have  been  evaluated  by  using  the  following  parameters 
(as  done  above):  dinc/dp  =  2.7  x  10  atm  =2.^xl0  dyne  cm  ; 

PTW  es  =  83*4x10"12  dyne"1  cm2;  e=7.0;  ^  =  14.7x10"  esucm; 
dfne/dT  =  5  x  10”3  K~  1,8/pcp  =  5.817  x  10  11cm3/erg.  The  results  for 
AV_  and  AV  are  reported  in  Table  V  for  various  collision  distances  a. 

*  U  J.  Q  Q 

Notice  that  usage  of  K„.  rather  than  K.,.  affects  AG  and  AS  ,  bpt  does 

0  Bj  Bj 

not  affect  AV_,  AH  and  hence  AVcas  reported  in  Table  V.  Comparing 

A  O 

the  results  of  Tables  IV  and  V  one  may  see  that  the  results,  especially 
because  of  the  AH  term,  differ  widely.  The  values  of  the  AV  ,  however, 

o  _  g  o 

seem  roughly  comparable  around  a  =  6-8x  10  cm.  They  also  are  of  the 

same  order  of  magnitude  as  the  experimental  AV.'s.  It  would  appear 

-8  b 

therefore  that  both  theories  suggest  6-8x10  cm  as  the  average  separa¬ 
tion  distance  of  the  quadrupoles  in  solution  of  1,2-DME.  Given  for  LiCfO 

-8 

the  separation  distance  of  the  charges  with  dipole  is  a  ,  =  3x  10  cm 

IQ  -  1  0 

(  |l/e  =  14.7x10  /4.8xl0  )  the  above  would  imply  a  "solvated" 


quadrupole.  The  situation  may  be  different  for  LiBF  in  2MTHF  given 

-8  ’ 
a^®»  2.5x10  cm  and  because  two  steps  of  association  to  quadrupoles 

have  been  detected  by  ultrasonic  relaxation  techniques  [  3  1 . 

For  the  case  observed  for  LiBF  .  in  2MTHF,  the  formation  of 

4 

quadrupoles  ought  to  be  written  according  to  the  scheme  [31 


2  LiB: 


Li  B  .  . .  Li  B 


(LiB). 


where  B  =  ligand  and  LiB  . . .  LiB  is  a  solvent  separated  dimer.  For  LiBF^ 
in  2MTH  the  formation  constant  for  the  second  step  of  the  above  scheme 
has  been  estimated  to  be  1^  =  0. 4.  The  calculation  of  the  volume 

change  associated  with  the  second  step  will  follow  the  same  lines  suggested 
by  Hemmes  [  14],  based  on  the  Kirkwood's  expression  [  4]  for  the  free 
energy  of  a  dipole  in  solution 

2 

CL"  3  2e  +1 
ri 

For  the  process  LiB  .  .  .  LiB^=5  (LiB^  +  S,  where  S  is  a  solvent  molecule,’ 
one  would  then  write 


AG°  = 


c  c.  c 

rMis  ^_s  Uos1  c  - 1 

"L  3  3  “  3  J  2c  +1  = 


-RT  in  K, 


where  the  subscript  "is",  "os"  and  "s"  stand  for  inner-sphere,  outer-sphere, 
and  solvent,  respectively 

^^2  _  rUis  ,  _h*os“l  d  ,  €-1  .  dc 


2,  T  dp 


-  r__LS  +  s  Mos  “I  d_  ,  e  -  1  .  de 
"13  _3  '  3  J  de  2c+l  ’  dp 


r.  r  r 
is  s  os 


giving 


r Oosi 
L  3  +  3'  3  J 


e  -  1 _ 3c  dine 

(2 e+ 1 )  (2e+l)(e-l)  dp 


=  -RT  inK 


2  (2c-H)(e-l) 


din  e 
dp 


-19- 


In  order  to  perform  the  calculation  of  AV,  _,  the  following  figures  have 

/  1  -5  - 1 

been  retained;  K^=0.4,  e=7,  dfne/dp  =  3x10  atm  .  This  gives 
AV2  ^,  =  0.16  cm  /mole.  The  same  sequence  as  above  will  now  be 
followed  for  the  calculation  of  AV2  g 

1Ar<o  2  2  2 

ac°  SAG2  rMis  ,  4s  uos]  d  (c-l)de 
A£>2  “  "  ST  "  L  3  +  3  “  r3  Jds  2c  +  1  dT 
rifl  rs  ros 


giving 


2  2 

^os  ~|  c - 1  _ 3s  dine 

3  ’  r3  J  (2c+l)  (2e+l)(e-l)  dT 

S  0  8 


*s2  =  RTlnK2  (g+iffr-u 


dine 

dT 


(34) 


Then 


-  AC”  +  TAS”  -  -  RT,nK2  +  RT2^  ^  . 


-RTinK2[l  - 


dine  *1 
din-T  J 


(35) 


The  calculations  of  AH2  with  the  same  parameters  as  above,  namely, 

K~  =  0.4,  c=  7,  dine/dT  =  5  x  10  gives  AH?  =  1 . 48  x  10  ^ergs/mole  = 

“  113 

1.48KJ/mole.  Then,  given  9/pc  «  5.8x10  cm  /erg,  it  gives 
.  3  P  , 

AH,0/ pc  =0.86  cm  /mole  and 
c  p 


AV,  _  =  AV_  AH-  =  0.  16-0. 86  =  -0.  70  cm^/mole 

it g  S  it p  1  pc  b 

P 


This  figure  appears  to  be  smaller  than  the  one  calculated  for  a  one- step 

-8 

process  for  a  =  6-8x10  cm.  The  same  conclusion  was  reached  by 
Hemmes  [  141  for  the  analogous  case  of  a  two-step  association  to  ion- 
pairs  starting  from  free  ions. 
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Table  I 


Dielectric  relaxation  parameters  eq,  e^  ,  f^,  and  specific 

conductivity  x  fr  LiBF^  in  1.2-DME  at  the  concentrations  investigated 
and  t  *=  25°C. 


0.092 

0.05 

0.035 

0.075 


c 

<x>2 

fl 

(GHz) 

f2 

(GHz) 

X 

(0  ^cm  *) 

2.6 

2.5 

35 

6.0  xl0~5 

2.7 

2.0 

35 

2.7 x  10"5 

2.6 

2.0 

40 

1.75x10“5 

2.7 

2.5 

40 

4.6  xlO“5 

Relaxation  parameters  A,  B  and  fr,  maximum  excess  sound  absorption  per 

wavelength  y  (  =^Afr  u  )  and  sound  velocity  u  for  LiBF.  in  1,2  -  DME 
m2  4 


at  the  concentrations  investigated  at  t'=25°C.  Calculated  y^'s  using 


Eqs. 


3  -1 

7  and  the  parameters  AV^  =  6.8  cm  /mole  and  K  =2. AM 


c 

(M) 

Ax  1017 

,  -1  2\ 
(cm  s  ) 

Bx  1017 

,  -1  2. 
(cm  s  ) 

fr 

(MHz) 

y  x  10"5 
m 

u  x  10"5 
(cm  s  *■) 

y  x  105 
m 

0.55 

51 

45 

428 

1.198 

421 

0.50 

55 

50 

19 

379 

1.197 

376 

0.A6 

53 

47 

■ 

382 

341 

0.415 

42 

52 

110 

277 

302 

0.38 

40 

51 

100 

238 

1.192 

271 

0.28 

36 

44 

95 

203 

1.190 

186 

x  n  x 

Values  of  b,  Kq,  K_,  K  ,  ,  Q,  and  as  functions  of  the  parameter  r  for 

o  o  *J  «J 

—18 

dipole-dipole  interaction,  where  y  =  14.7  x  10  esucm,  e  =  7.0,  T  =  298K. 


r  x  10 


0.4343 

0.5637 

0.7503 

1.029 

1.465 

2.187 

3.474 

6.002 


K 

o 

M-1 

*F 

m"1 

4.357 

6.73 

3.356 

5.90 

2.522 

5.34 

1.838 

5.14 

1.291 

5.59 

0.865 

7.71 

0.545 

17.6 

0.315 

127.4 

.590 
.454 
0.341 
0.249 
0.175 
0.11? 
0.0737 


0.723 

0.696 

0.756 

1.044 

2.377 

17.27 


0.3234 

1.1933 

2.1828 

3.5448 

6.3303 

20.597 


V 

M**1 

* 

Vj 

M-1 

0.61 

0.082 

2.26 

0.306 

4.13 

0.559 

6.71 

0.908 

11.98 

1.62 

38.97 

5.27 

Notice  that  following  Eq.  20  K_  =K  bQ,  Q  has  been  calculated  with  the  series 

01  o 


expansion  up  to  n=10.  Also  q  =  10.40x10  cm  in  this  case. 


Table  IV 


Calculated  AV^,  AHQ  and  AVg  for  a  dimerization  of  dipoles  in  1,2-DME  at  T=  298.2K 
according  to  Eqs. 14,  25  and  27. 


ax  10 


RTb  d*nE 
-  dp 

cmJ/mole 

RTB 

cm J /mole 

avt 

cm-Vmole 

a 

AHo 

KJ/mole 

(0/pCp)AHo 

cm^/mole 

b 

AVg 

o 

cnr/mole 

0.28 

2.07 

-1.79 

-2.68 

-1.56 

-0.23 

0.37 

2.07 

-1.70 

-3.48 

-2.02 

0.32 

0.49 

2.07 

-1.58 

-4.63 

-2.69 

1.  11 

0.67 

2.07 

-1.40 

-6.35 

-3.70 

2.30 

0.96 

2.07 

-1.11 

-9.05 

-5.26 

4.  15 

1.43 

2.07 

-0.64 

-13.5 

-7.86 

7.22 

2.26 

2.07 

0.  19 

-21.5 

-12.5 

12.7 

3.91 

2.07 

1.84 

-37.1 

-21.6 

.23.4 

'7.64 

2.07 

5.57 

-72.4 

-42.1 

47.7 

0 


Table  V 


Evaluated  AV^.  AH0  and  AVS  for  a  dimerization  of  dipoles  in  1,2-DME  atT=  298.2K 
according  to  Eqs.  29  and  31. 


ax  10** 
(cm) 

avt 

cm"* /mole 

a 

AH0 

ergs/mole 

(0/ pCp)  AH0 
cm^/mole 

b 

avs 

cm  /mole 

10 

4.32 

-6.00  x  1011 

-34.90 

39.2 

9 

0.09 

-2.02  x  1011 

-11.77 

11.9 

8 

-0.52 

-1.45  x  1011 

-  8.45 

7.93 

7 

-0.65 

-1.33  x  1011 

-  7.74 

7.09 

6 

-0.45 

-1.52  x  1011 

-  8.84 

8.39 

5 

+0.74 

-2.64  x  1011 

-15.35 

'  16.  1 
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Figure  Captions 

Fig.  1  a.  Log10ft  vs.log10  c  for  LiBF4in  l,2-DMEat25°C. 

b.  ilgfciyTvs.  c(l  --flfar  LiBF4in  1.2-DME  at  25°C. 

*'o 

Fig  2  Representative  plot  of  e  '  and  e"  vs.  f  for  LiBF^  in  1,2  DME  at 
t=  25°C. 

Fie  3.  Representative  plot  of  the  quantity  Ct/f2  vs.  the  frequency  f  for 
LiBF^in  1,2-DME  at  25°C. 

Fig.  4  Bottcher  plot  for  LiBF4  in  1,2-DME  at  t  =  25°C. 

• 

Fig.  5  a.  vs.  c  for  L-iBF^  in  1,2-DME,  t  =  2  5°C. 

b.  (TTC2/2gs^mRT}  vs.  c  for  LiBF4  in  1 , 2- DME  at  25°C. 

c.  Calculated  (solid  line)  and  experimental  values  of  pm  vs.  c 
for  LiBF4  in  1,2-DME  at  t  =  25°C. 

Fig.  6  Plot  of  the  functions  K^,,  K^,  and  Kgj  vs.  tbe  dipole-dipole 

separation  distance  in  a  medium  of  permittivity  e  =  7  for  dipoles 
of  moment  (j  =  14.  7  x  10“  esu  cm  at  T  =  298 ,2K. 


•Mi 


0.02  0.04  0.06  0.08 
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